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ABSTRACT 

In the context of the unified model of Seyfert galaxies, we use observations from 
the literature and a radiative transfer model to investigate the near-IR to mm emission 
produced by the presumed torus in the Circinus galaxy, from 2/im to 1.3mm. From 
the infrared SED modelling, we find that the total luminosity (L//?) in this wavelength 
range, consists of similar contributions from the torus and starburst with a ratio of 
nuclear luminosity to starburst luminosity (Lnug/^Sb) ~ 0.8. 

By using a similar torus model to that of NGC1068, but without the conical dust, 
we find an upper limit to the outer torus radius of ^12pc with a best fit of ^2pc. The 
upper limit torus size estimated from the radiative transfer modelling is consistent 
with the 16pc torus radius estimated from near-IR imaging polarimetry of Circinus. 

Key words: galaxies: individual: Circinus - galaxies: nuclei - galaxies: active - in- 
frared: galaxies - radiative transfer 



1 INTRODUCTION 

Observational evidence for obscuring material in the centres 
of Active Galactic Nuclei (AGN) is now abundant. Obscu- 
ration is observed in the form of molecules and dust, present 
in the centres of these galaxies on scales that vary from a 
fraction of a parsec to tens and even hundreds of parsecs 
(Malkan et al. 1998; Gallimore; Siebenmorgen et al. 1997; 
Greenhill et al. 1996). It is believed that the obscuring ma- 
terial takes the form of a dusty "torus" which blocks and 
absorbs part of the radiation from the active nucleus. These 
dusty tori play an essential role in unified theories of Seyfert 
galaxies and strongly support the unification models. 

The unified model for Seyfert galaxies proposes that all 
types of Seyfert galaxy are fundamentally the same, but, 
the presence of the dusty torus obscures the nucleus, includ- 
ing the broad line emission region in many systems. In this 
picture the classification of Seyfert 1 or 2 depends on the 
inclination angle of the torus to our line of sight (e.g. An- 
tonucci and Miller 1985; Antonucci 1993). Strong support for 
this model has come from X-ray observations, showing that 
Seyfert 2 galaxies are Compton thick, i.e. the nuclear radia- 
tion is absorbed by matter with column densities >10^*cm~^ 
(Lawrence and Elvis 1982; Maiolino et al. 1998a). 

Direct imaging of the presumed torus is technically de- 
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manding due to its (predicted) small size (e.g. Efstathiou, 
Hough and Young 1995, hereafter EHY95; Granato, Danese 
and Franceschini 1997; Alexander et al. 1999). The most 
convincing direct evidence for a Seyfert torus comes from 
observations of the HON molecule (Jackson et al. 1993), 
which traces dense molecular gas, and near-IR imaging 
polarimetry (Young et al. 1996), both of NGC1068. These 
images show a large nuclear structure of ~200pc in ex- 
tent, approximately perpendicular to the [Olll]5007 emis- 
sion cone and radio lobe emission. Indirect support for 
tori comes from ground-based narrow band imaging (e.g. 
Wilson and Tsvetanov 1994), high resolution HST imag- 
ing (e.g. Capetti et al. 1997; Falcke, Wilson and Simpson 
1998) and imaging-polarimetry (e.g. Ruiz et al. 2000; Lums- 
den et al. 1999; Tadhunter et al. 1999; Packham et al. 1997) 
which have shown a number of objects with nuclear cone- 
like structures, presumably collimated by the torus. 

The origin of the infrared continuum in a number of 
Seyfert 2 galaxies is clearly thermal as indicated by sili- 
cate features such as the broad absorption at 10/^m (Roche 
et al. 1991) and the various emission features from poly- 
cyclic aromatic hydrocarbons (PAHs), e.g. 3.3, 6.2 and 
11.3/im, which are normally associated with the presence 
of hot stars in a starburst environment. Dust grains, dis- 
tributed in a toroidal geometry surrounding the central 
black hole and accretion disc, are believed to be respon- 
sible for the absorption of radiation from the central source 
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and re-radiation of energy into longer wavelengths giving 
rise to the observed "IR bump" longwards of Ifim and peak- 
ing at mid-IR wavelengths (Rowan-Robinson and Crawford 
1989). The resultant radiation spectrum is dependent on the 
distribution of dust grain temperatures (< 2000K) and op- 
tical depth. Two additional components can add to the IR 
spectra of active galaxies (Rowan-Robinson and Crawford 
1989): galactic diffuse emission from large dust grains heated 
by the interstellar radiation field in the galactic disc (cirrus, 
Tiarge-grains < 40K ) and a starburst component peaking 
at about 60/im. 

Thus, a radiative transfer model which takes into ac- 
count dust grain features within a true toroidal configuration 
and starburst environment is needed to correctly reproduce 
the near-IR to mm SED of the Circinus galaxy. 

Here, we present a model for the IR continuum emission 
of the Circinus galaxy which successfully fits the observa- 
tions from 2fj,m to mm wavelengths. We constrain the size 
and flux of the dusty AGN torus and estimate the contribu- 
tion of the starburst emission to total IR emission. Section 2 
presents the data compilation and sections 3 and 4 present 
the model and the fitting procedure. Section 5 discusses the 
model and its implications and conclusions are summarised 
in section 6. 



2 THE DATA 

The continuum and PAHs data points have been compiled 
from a variety of sources in the literature and are presented 
in Table 1. The data sources will be described below. There 
is a large range of sizes in the observational apertures, from 
arcsec to arcminute beamsizes, equivalent to physical sizes 
of <5pc and >200pc (for Circinus, 1 arcsec corresponds to 
20pc); whenever possible, data from the smallest available 
aperture were used to minimise the contribution from galac- 
tic disk emission. 

In the near-IR, below 2/im, the stellar contribution is 
quite significant, and to study the nuclear emission, the 
starlight contribution has to be subtracted from the inte- 
grated flux. The data point at 2.2^m shown in Table 1 
corresponds to the deconvolved non-stellar measurement 
as presented by Maiolino et al. (1998b) and derived from 
a O.lSarcsec aperture. This measure will set a limit on the 
size of the nuclear non-stellar source (see section 2.1). 

In the mid-IR, between 3 and 20^m, the continuum 
data points at 3.8 and 4.8/im, are provided by Siebenmor- 
gen et al. (1997) who carried out speckle interferometry in 
the L'(3.8/im) and M(4.8/im) bands with a 0.3arcsec aper- 
ture. The speckle observations are insensitive to larger scale 
components, e.g. starlight, and predict more directly the nu- 
clear non-stellar contribution (see section 4). Other contin- 
uum data points over this wavelength range are from the 
ISO-SWS spectrum presented in Sturm et al. (2000). 

Line information for the PAH features at 3.28, 6.2, 7.7, 
8.6 and ll.S^m were directly measured from the ISO-SWS 
spectrum in Sturm et al. (2000). The reported values in Ta- 
ble 1 correspond to the peak of these emission lines. ISO 
apertures are large enough (~24arcsec'^) to contain much of 
the circumnuclear star forming regions as well as the nuclear 
emission. The ISO spectrum shows that the region between 
5 and 12^m is dominated by PAH features, characteristic 
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nuclear measurement 

PAH feature peak 
" 5a upper limit 

of star formation activity (see Fig 1), and it is remarkably 
similar to those of pure starburst galaxies also observed by 
ISOPHOT (Siebenmorgen, Kriigel and Zota 1999). 

The spectrum corresponding to the Si absorption fea- 
ture centered at 9.7/^m is shown in more detail in Fig 
1. It was provided by P. Roche and published by Roche 
et al. (1991). This spectrum was taken with the UCL spec- 
trometer and used an aperture of 4.3arcsec. For compar- 
ison also shown in Fig 1 is the mid-IR ISO-SWS spec- 
trum of Circinus taken with the much larger aperture of 
14x20arcsec^ (Sturm et al. 2000). 

The N(10.3/im) band continuum data point, is also pro- 
vided by Siebenmorgen et al. (1997) and measured through 
a 1.3arcsec slit. It is interesting to note that a larger aper- 
ture measurement at 10.3/^m by Moorwood and Glass (1984) 
with a 5arcsec aperture is very similar to the smaller aper- 
ture data of Siebenmorgen et al. (1997). This is a good 
indication of the compact nature of the central source at 
this wavelength. The continuum data points at 20, 30 and 
40^m were taken from the ISO-SWS spectrum (Sturm 
et al. 2000). The ISO SWS spectrum at 25/im agrees within 
20% of the IRAS flux. 

The far-IR continuum observations were collected from 
the IRAS large aperture observations (OOarcsec) at 12, 25, 
60 and 100/xm as reported by Moshir et al. (1992). The con- 
tinuum data point at 150/im wEis taken from the analysis 
by Ghosh et al. (1992) which corresponds to an aperture 
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of 40arcsec. This measurement is only a 5a upper limit as 
Circinus was not detected at this wavelength. 

At millimeter wavelengths, Circinus has been detected 
at 1.3mm with an aperture of 23arcsec as reported by 
Siebenmorgen et al. (1997). Maps of the millimeter contin- 
uum emission show Circinus as an unresolved source. 

Clearly, the large aperture observations are contami- 
nated with light from various emission sources particularly 
from the extended galactic disc; thus, special care is needed 
when fitting the large aperture data points to the individ- 
ual emission components in the radiative transfer model (see 
next section). 

2.1 Constraints from Observational Data 

In the K band (2.2/im), the source is clearly unresolved in 
a O.Sarcsec aperture (6pc on source). By deconvolving the 
nuclear K-band surface brightness radial distribution with 
a stellar profile and a nuclear PSF, Maiolino et al. (1998b) 
set a constraint on the size of the nuclear non-stellar source 
of <1.5pc. Similarly, at 3.8 and 4.8^m the source is also ob- 
served to be unresolved at 0.7arcsec resolution and is possi- 
bly <0.3arcsec in size, putting upper limits of 6pc (Sieben- 
morgen et al. 1997). Clearly, at wavelengths <5/im most of 
the emission comes from a central source with a size <6pc. 
This size scale almost certainly corresponds to the inner ra- 
dius of the proposed dusty torus, as the near-IR emission is 
most likely to result from hot dust grains (T<1300K) heated 
by the intense radiation from the non-stellar central source. 
Thus, in our model, the near-IR data points at <5jj,m are 
assumed to be totally dominated by the emission from the 
torus component. 

Observations at lO^m reveal that the central source size 
is possibly resolved with a size of 26pc whilst at 20^m it is 
unresolved and <30pc in size (Siebenmorgen et al. 1997). 
As mentioned in the previous section, the emission detected 
with the larger aperture of Moorwood and Glass (1984) at 
lO^m does not seem very different from that of smaller aper- 
ture by Siebenmorgen et al. (1997). At wavelengths larger 
than 20fim, the IR emission most likely results from a more 
extended region, and there is possibly a significant contri- 
bution from the inner circumnuclear starburst (~ 40pc in 
radius) . 

The spatial resolution analysis at 60, 100 and 150/im by 
Ghosh et al. (1992) reveals a deconvolved size of Circinus of 
40arcsec at 50/im. However, at these longer wavelengths, the 
resolution is too poor to put useful constraints on the size 
of the emitting source. 

The spectrum corresponding to the Si feature at 9.7/im 
(Fig 1) was taken with an on-source equivalent aperture of 
~80pc (Roche et al. 1991). The detection of a number of CO 
stellar absorption features in a 4.4x6.6arcsec^ aperture gives 
strong evidence for circumnuclear star-formation on a re- 
gion few tens of parsecs from the nucleus (Oliva et al. 1994, 
Maiolino et al. 1998b), thus, the Si feature cannot be solely 
modelled by dust self-absorption from the torus component, 
but, a starburst contribution will have to be included when 
fitting this feature with the SED model. 

The long wavelength data points and PAH features 
(IRAS and ISO data), which correspond to large aperture 
observations (see Table 1), are hence dominated by the 
larger scale starburst emission observed in the circumnuclear 



ring at ~200pc from the nucleus. Thus, the data points cor- 
responding to IRAS and ISO data are fitted with a starburst 
component. 

3 THE GENERAL INFRARED MODEL 

The model used is the radiative transfer model of Efstathiou 
et al. as described in detail in Efstathiou and Rowan- 
Robinson (1995, hereafter ER95), EHY95, and Efstathiou, 
Rowan-Robinson and Siebenmorgen (2000, hereafter ERS). 
This radiative transfer model has been successful in fitting 
the IR continuum of NGC1068 (EHY95) and Centaurus A 
(Alexander et al. 1999). The proposed model for the SED 
of the Circinus galaxy is a combination of the following two 
components described below. 

3.1 Torus 

As concluded in ER95, the most likely geometry for the ob- 
scuring material in the center of AGNs, is that of a "tapered 
disc". Other proposed geometries for the torus are "flared 
discs" (Granato and Danese 1994; Efstathiou and Rowan- 
Robinson 1990 (hereafter ER90)) or cylinders (Pier and Kro- 
lik 1992; 1993; Taniguchi and Murayama 1998). However, 
the shape of the IR continuum spectrum of type 1 and type 
2 objects, the appearance of the lO^m silicate features and 
the statistics of the two types of Seyfert galaxies, give strong 
support for a tapered disc geometry (see ER95 for detailed 
discussion). Indeed, there is further confirmation of the va- 
lidity of the tapered disc geometry to realistically represent 
the nuclear emission of tori in AGNs. This comes from the 
fitting of IR data to a sample of Seyfert galaxies with the 
EHY95 models of Figure 3. They have been successfully used 
to fit nuclear (non-stellar) near-IR observations for a sam- 
ple of Syls and Sy2s (Alonso-Herrero et al. 2001) showing 
that the near-IR nuclear emission from Seyfert galaxies can 
be fitted with the tapered disc geometry tori of EHY95. We 
therefore assume this geometry for the Circinus torus. As 
demonstrated in ER90, the orientation of the plane of the 
system to the line-of-sight is a very important parameter 
when fitting the observed spectra, since only a small differ- 
ence in this inclination parameter, dv, produces a significant 
variation in the emergent spectrum. Another important pa- 
rameter that determines the shape of the IR continuum, is 
the half-opening angle of the torus, O. In general, these pa- 
rameters are not known for most objects, but fortunately, 
we have a good idea of their values from spectropolarime- 
try (Alexander et al. 2000) and IR-imaging polarimetry of 
Circinus (Ruiz et al. 2000) which gives Sv~40 degrees (mea- 
sured from the axis of the torus to the line of sight) and 
0~45 degrees. Although the emission line cones observed 
in Circinus and many other AGNs offer a measure of the 
opening angle of the torus, it is common to find that the 
conical emission "opens-up" as light travels and scatters in 
the circumnuclear medium. Thus, we should keep in mind 
that O can be smaller than that derived from emission line 
cones (see section 4). 

Other parameters that will affect the shape of the emer- 
gent SED spectrum are: the ratio of the inner and outer radii 
(ri/r2), the ratio of the height to the outer radius {h/r2), 
the equatorial optical depth C^'Ajjv), the dust sublimation 
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Wavelength (microns) 



Figure 1. Circinus Mid IR spectrum. Solid line corresponds to the Circinus ISO-SWS (Sturm et al. 2000) spectrum and square open 
points correspond to the Roche etal data. The ISO spectrum is dominated by emission features from PAHs as the ISO-SWS spectrum 
was taken with a large aperture (14x20arcsec^ ). Note the lack of these emission features in the Roche etal data (aperture ~ 4.2arcsec). 



temperature (Ti) and the radial dependence of the density 
distribution r"''. To take into account the effects of dust 
in the nuclear regions of Circinus, we assume that there is 
an additional extinction (i.e. in addition to that produced 
by the torus). Marconi et al. (1994) found Ay > 20mag to- 
wards the nucleus. This limit follows from the fact that we 
can see the nucleus at 1.25/im but are unable to do so at 
7000A, thus, they only set a lower limit for the nuclear ex- 
tinction. However, we find that Av-=25 mag is required to 
successfully fit the near-IR data. This value is also close to 
that found by Moorwood and Glass (1984) from the depth of 
silicate feature at 9.7/im and is also consistent with that es- 
timated by Maiolino et al. (1998b). Based on H and K band 
observations, they set a minimum extinction Av > 12mag 
toward the nucleus of Circinus. From their images, it is clear 
that this obscuration is caused by the large amounts of dust 
as seen in their H-K colour map and it is due to the presence 
of a nuclear gas bar. Their estimated Ay however is only a 
lower limit and we expect a larger value since the extinction 
is most likely to be due to an inhomogeneous distribution of 
dust, gas and stars. The parameters assumed for the torus 
model are listed in Table 2. 



3.2 Starburst 

The radiative transfer models of starburst galaxies are de- 
scribed in detail in ERS. The models consider a starburst 
galaxy as an ensemble of optically thick giant molecular 
clouds centrally illuminated by recently formed stars. The 
stellar population is determined according to the Bruzual 
and Chariot (1993) stellar population synthesis models. The 
grain model used is that of Siebenmorgen and Kriigel (1992) 
and the emission of transiently heated dust grains is calcu- 
lated according to their method. In the models, the effects 
of multiple scattering is taken into account. The models re- 
late the observed properties of a galaxy to its age and star 
formation history. These models have successfully matched 
the observational characteristics of M82 and other starburst 
galaxies (ERS). Various populations of dust particles are 
considered: 

(i) Mathis, Rumpl and Nordsieck (1977). These are large 
particles with a size distribution of classical grains 

(ii) small graphite grains, to account for the extinction bump 
at around 2 175 A (Draine 1989) 

(iii) PAH molecules, to explain the near-IR and mid-IR 
emission bands 

To take into account the effects of dust in the starburst 
regions of Circinus, we extinct the starburst emission by 
Av=5mag, as measured from narrow emission line ratios 
(Oliva et al. 1994). This is lower than the additional extinc- 
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Best fit torus parameters ^1/^2 0.05 

h/r2 0.5 

{°) 50 inner radius (ri, pc) 0.12 

©I (°) 30 outer radius (r2, pc) 2.4 

1000 height (h, pc) 1.2 



system inclination|j, 0v, 
cone opening half-angle, 
"''Auv (mag) 

'^lAv (mag) 200 
'""Ay (mag) 830 
TUKJQ 1000 

^ P ' 

additional AyMmag) 25 



Best fit starburst parameters 

Starburst e-folding time 10 Myr 

Starburst age 26 Myr 

TV, initial GMC optical depth 50 

Ay in ISM (mag) 5 

Ay through Galactic disk (mag) 1.5 



measured from the plane of the torus to the line of sight 
^ '^'^Ay corrected for the inclination of the torus to the line of sight 
dust sublimation temperature at torus inner radius 
value determined as discussed in EHY95 
due to circumnuclear bar 

tion assumed to the nucleus (section 3.1), as the starburst 
is expected to be outside the circumnuclear molecular bar. 



4 MODEL FITTING 

As described in section 3, the data points will be fitted to 
the SED model components according to the observational 
constraints described in section 2.1. The critical wavelength 
range is likely to be where the torus and starburst compo- 
nents are equally important, that is in the 10-30/im region. 
Here, the strong absorption feature is indicative of the pres- 
ence of dust in the form of silicates which are likely to be 
present in a dusty environment such as the torus, however, 
the observational beam for this data is large enough to cover 
a region which also includes star formation. Thus, contribu- 
tions from both the torus and starburst components are im- 
portant in this wavelength range. On the other hand, given 
that the near-IR points are to be modelled by the torus 
alone (see section 2.1), this sets a strong constraint on the 
number of possible models to be tested for the torus. Simi- 
larly, in the long wavelength range and large aperture data 
points which include the PAHs features, the starburst com- 
ponent dominates, and only a number of starburst models 
can be applied. Thus, the overall fitting procedure is much 
better constrained and the number of free parameters re- 
duced once the short and long wavelength ranges are con- 
sidered. 

The torus input parameters that reproduce the near- 
IR points are listed in Table 2. The best fits were obtained 
with a cone half-opening angle of 30 degrees, the same as 
the best fit obtained for the NGC1068 torus model (EHY95). 
One important difference with the NGC1068 model is that 
there is no requirement to include a contribution from hot 



dust in the torus cone. Indeed, if this is included, then it is 
very difficult to obtain a fit to the near-IR data points. The 
starburst parameters which reproduce the PAHs features 
and large scale emission successfully are also listed in Table 
2. Table 3 presents the derived geometrical characteristics 
for the best fit torus in Circinus. In Figure 2 we show the 
best model for the SED of the Circinus galaxy which is the 
addition of the torus and starburst models. Note that the 
solid line (representing the total emission from the galaxy) 
exceeds the mid-infrared data of Roche et al. (1991). This is 
not surprising as the starburst which contributes about half 
of the emission is more extended that the 4arcsec aperture 
of Roche et al. (1991). Thus, the final SED for the lO^m 
absorption feature should not fit the 9.7/im data of Roche 
et al. 1991 (Fig 1). 

Unfortunately, there is a serious lack of small aperture 
(nuclear) data in the range of 30-60/im. As a consequence, 
we are unable to rule out tori up to ~12pc although in- 
creasing the torus outer radius further would produce far 
too much 30-60/^m flux. Future observations in this critical 
wavelength range, will put extra constraints on the torus 
model, and the torus size would be determined more accu- 
rately. 

It is interesting to note in Figure 2 the overall larger 
contribution to the total IR emission of the starburst com- 
pared to the torus component. This is indicated by the ratio 
of the nuclear luminosity -torus component- to starburst 
luminosity -starburst component- {Lmuc / Lsb) ~ 0.8 (see 
next section). 



5 DISCUSSION 

Our SED model shows that, overall, the starburst contri- 
bution to the total IR luminosity, Lm slightly dominates 
over the torus contribution. A direct result of this, is the 
relatively small size of the torus in Circinus, with a best fit 
of ~2pc for the outer radius, comparable to the torus outer 
radius of 1.8pc modelled for Cen A (Alexander et al. 1999) 
but significantly smaller than the ~90pc torus for NGC1068 
(Young et al. 1996 and EHY95). This result is smaller than 
that calculated from the near-IR imaging polarimetry mod- 
elling of Ruiz et al. (2000) which suggests that the torus is 
approximately 16pc in radius, although consistent with our 
upper limit of ~12pc. 

So far, by using the radiative transfer models of Efs- 
tathiou et al. we have been able to model the size of the 
torus outer radius for three objects: NGC1068 (EHY95), 
CenA (Alexander et al. 1999) and Circinus (this work). Al- 
though the statistics are small, it is interesting to note that 
all three AGN can be modelled with tori of very similar pa- 
rameters. This suggests that whatever the mechanism that 
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log (h/fim) 

Figure 2. IR SED for the Circinus galaxy. The SED is fitted with a starburst and a torus model. Dark circles are the data points 
from Table 1 and Figure 1. Units for i/Sj^ are 10~^^W/m^. Dotted line corresponds to the torus component, dashed line is the starburst 
component and the solid line represents the total fit. The model parameters assumed for the torus and starburst are given in Table 2. 



forms these tori in an AGN-SB environment, it is such that 
the size of the tori scales with the luminosity of the cen- 
tral source but the rest of the torus parameters (Table 2) 
are scale invariant. For example, the hard X-ray contin- 
uum, which is considered to be a good indicator of the cen- 
tral source strength, indicates that the Circinus galaxy is 
a relatively low luminosity AGN (Matt et al. 1996), while 
NGC1068 is ~ 100 times more powerful (Turner et al. 1997), 
with a torus radius ~100pc. The predicted torus radii scales 
as L^^^ , reflecting the fact that both objects can be fitted 
with very similar radiative transfer models. 

Of great debate regarding the nature of the Circi- 
nus galaxy is the ratio of the AGN to starburst emission 
{Lnuc/Lsb), that is, the relative contribution of the torus 
and starburst components to the total Lm. Previous predic- 
tions from IR modelling differ widely. As already mentioned, 
Rowan-Robinson and Crawford (1989) found a 20% contri- 
bution of the nuclear component to the total i/ji(io-iooMm) 
and Lnuc/Lsb ~0.3 while this work suggests a Lnuc/Lsb 
of 0.8. On the other hand, from the Bra luminosity, Moor- 
wood et al. (1996) estimate that this ratio is perhaps as high 



as 9. Our result is in very good agreement with that sug- 
gested by Maiolino et al. (1998b) who estimate from K band 
observations a Lnuc/Lsb ratio of 0.9 for the region <200pc. 
Similarly, we agree with the ratio derived from the inferred 
ionising continuum required to produce the highly ionised IR 
lines suggesting a ratio closer to 1 (Moorwood et al. 1996). 
We believe that the differences are caused by a number of 
pitfalls. For example, Rowan-Robinson and Crawford (1989) 
modelled the far-IR (lO-lOO/im) IRAS data with a three 
component model, a starburst, a Seyfert and a disc cirrus- 
like component. This model assumed a spherically symmet- 
ric dust distribution with a power-law density distribution, 
n(r) ~ . In this model, radiation from the polar regions 
of the torus is absorbed by the spherical dust cloud. Thus, 
this model will naturally under-estimate the nuclear com- 
ponent, as their low Lnuc/Lsb ratio shows. However, there 
is now ample observational evidence that favours the exis- 
tence of axially symmetric distributions of dust (Efstathiou 
and Rowan-Robinson 1990, ER95). On the other hand the 
large Lnuc/Lsb ratio derived by Moorwood et al. (1996) 
suggests a strongly dominant AGN over the starburst. This 
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was derived from the comparison of near-IR recombination 
lines in Circinus and M82; from a Br7 map (Moorwood and 
Oliva 1994), they assume that 50% of the recombination 
lines are produced in the nucleus of Circinus. However, this 
is likely to be an over-estimation of the strength of Br7 in 
the nucleus since,, more recently, a more accurate measure- 
ment of the Circinus nuclear recombination lines (Maiolino 
et al. 1998b) has shown that only ^^10% of the total recom- 
bination lines flux is produced in the nucleus of Circinus. 
Thus, assuming the standard relation between Bra flux and 
total IR luminosity in starburst galaxies, we can determine 
the contribution of the starburst to L/_r. We find Lsb 
8.5x10®Lq, so Lnuc / Lsb ~ 1-2 for an AGN luminosity of 
10^° Lq. This is in good agreement with the ratio we find 
from the SED modelling. 

Radio observations suggest that at these wavelengths, 
the starburst emission is the dominant component (Forbes 
and Norris, 1998) and Circinus would certainly be classified 
as a starburst galaxy based purely on radio emission. Indeed, 
it is interesting to note that the FIR-radio correlation estab- 
lished by Helou, Soifer and Rowan-Robinson (1985) which 
is followed by normal spiral and starburst galaxies, is closely 
followed by the Circinus galaxy, confirming Forbes and Mor- 
ris result. However, it is almost certain that this correlation 
is the result of large scale star-formation activity unrelated 
to the Seyfert activity which generates both the synchrotron 
radio emission and the thermal FIR emission. Given that the 
Circinus galaxy has strong star-formation activity on small 
and large scales, it is not surprising that this galaxy fol- 
lows the FIR-radio correlation. Radio observations on arc- 
sec resolution can not easily distinguish between Seyfert and 
starburst activity. The radio emission from both starburst 
regions and Seyfert nuclei have similar spectra and morphol- 
ogy, and the steep spectrum of even the compact cores seen 
in VLA maps of Seyferts suggests that they might contain a 
significant nuclear starburst component (Norris et al. 1992). 
On the other hand, long-baseline radio interferometry is sen- 
sitive only to compact, high brightness objects and is an 
ideal tool to discriminate over the two spatial scales. Unfor- 
tunately no radio interferometry is available for the Circinus 
galaxy and we can not clearly determine the nature of the 
central source from present radio data. 

Nevertheless, there is a clear large variation in the 
predicted Lnuc / Lsb ratios for Circinus, from starburst- 
dominated at FIR wavelengths, to AGN-dominated in 
the near and mid IR. Indeed, NGC1068, the prototypical 
Seyfert 2 galaxy shows a similar behaviour. The mid-IR(5- 
16/im) spectrum of NGC1068 is ^85-95% dominated by the 
AGN, and Lnuc /Lsb ~ 3. In the FIR however, where 
extended emission dominates, the nucleus does not con- 
tribute more than 25% to the total IR flux at 450/im (Le 
Floc'h et al. 2001). Clearly, determination of Lnuc / Lsb is 
strongly dependant on the spectral region where it is calcu- 
lated. 

Our SED modelling has the advantages over other pre- 
dictions of having a more complete model for the SED, in- 
cluding the emission of PAH features and a more realistic 
torus model. Additionally we have used a large set of re- 
cently published data, from NIR to mm wavelengths. More 
importantly, we pay special attention to the fit of data points 
from the various aperture sizes which give important con- 
straints to the size of the emitting source. The undoubted 



composite nature of the Circinus galaxy is shown in our IR 
SED model which predicts a galaxy dominated by the AGN 
torus at near and mid-IR wavelengths while at far-IR to 
mm wavelengths the starburst activity is the dominant one. 



6 CONCLUSIONS 

We can summarise our results as follows: 

• The case of the Circinus galaxy provides strong support 
for the dusty torus model of active galaxies which further 
support the unification scheme, that is, thick dusty tori with 
a tapered disc geometry surrounding the central engine. 

• Our best fit for the outer radius of the torus in Circinus 
is ~2pc, but a radius as large as 12pc can not be ruled out. 

• We confirm the composite nature of the Circinus galaxy, 
which is dominated by the AGN torus emission in the near- 
IR but largely dominated by starburst phenomena in the 
far-IR and mm wavelengths. The Ljr is slightly dominated 
by the starburst emission and the ratio Lnuc / Lsb ~ 0.8. 
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